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Abstract 


The Finite Difference Time Domain (FDTD) technique has been 
applied to a wide variety of electromagnetic analysis problems, 
including shielding and scattering. However, the method has not 
been extensively applied to antennas. In this short paper 
calculations of self and mutual admittances between wire antennas 
are made using FDTD and compared with results obtained using the 
Method of Moments. The agreement is quite good, indicating the 
possibilities for FDTD application to antenna impedance and 
coupling. 


I Introduction 

The Finite Difference Time Domain (FDTD) technique has had only 
limited application to antennas. This is somewhat surprising, 
since the geometrical and material generality of the method 
suggests that it might have significant application to antenna 
analysis, especially in situations where other structures, 
especially electromagnetically penetrable ones, are nearby. This 
is due to the relative ease with which the FDTD method 
accommodates modeling of volumetric electromagnetic interactions 
with materials as compared to the Method of Moments. 

Earlier work [1] has shown that the FDTD method could compute the 
self impedance of a wire antenna in three dimensions, however, 
the approach used in [1], plane wave incidence, did not lend 
itself to mutual coupling calculations. Accurate self-admittance 
FDTD results for two-dimensional antenna geometries were 
presented in [2]. In this paper we demonstrate both self and 
mutual admittance FDTD calculations for three dimensional wire 
antennas . 
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II Approach 


The test problem geometry is shown in Figure 1. Two wire dipoles 
of length 57 and 43 cm are parallel and separated by 10.5 cm. 

Both are center fed, and are symmetrically positioned. The goal 
is to determine the self admittance of the driven dipole and the 
mutual admittance between the two dipoles. 

The FDTD computations were made using a three dimensional 
computer code based on the Yee [3] cell, with second order Mur 
[4] absorbing boundaries. The problem space was chosen as 
61x51x80 cells, with the cell dimensions Ax=Ay=0.5 cm, Az=1.0 cm. 
Making the two transverse dimensions smaller results in a greater 
length to diameter ratio, so that a thin wire Moment Method code 
may be used to provide comparison results over a wider band of 
freguencies. Thinner wires may be modeled in FDTD using sub-cell 
methods [5,6] . 

For the FDTD calculations the longer dipole is fed at the center 
with a Gaussian pulse of 100 volts maximum amplitude that reached 
its 1/e amplitude in 16 time steps. The time steps were 11.11 
picoseconds, the Courant stability limit for the cell size 
chosen. 

During the progress of the FDTD calculations the currents at the 
center of each dipole were saved for each time step. They were 
computed by evaluating the line integral of the magnetic field 
around the dipoles at the center. Along with the applied 
Gaussian voltage pulse the currents were Fourier transformed to 
the frequency domain. Then, based on the admittance parameter 
equations 


h = v i 


U 21 


V Y 
v 2 *22 


and taking V 1 to be the driven dipole voltage and V 2 zero, we 
eas'ily obtain the self admittance of dipole 1 and the mutual 
admittance (since Y 12 = Y ?1 ) between the dipoles by dividing the 
appropriate complex Fourier transforms of V,, I H and I 2 . 

The Moment Method results were obtained using the Electromagnetic 
Surface Patch Version 4 [7] computer code. The wire radius for 
the Moment Method calculations was taken as 0.281 cm, providing 
the same cross section area as the 0.5 cm square FDTD cells. 

While the FDTD calculations should be valid up to approximately 3 
GHz based on having 10 FDTD cells per wavelength, the thin wire 
approximation for the Moment Method code becomes questionable at 
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approximately 1 GHz and this was taken as the upper frequency 
limit for comparison of results. 


Ill Results 

Figure 2 shows the Gaussian pulse voltage applied to the 1 cell 
gap at the center of the longer, driven dipole. Figures 3 and 4 
show the current flowing in the center cell of the driven and 
passive dipole respectively. All are plotted on the same time 
scale, corresponding to 8,192 time steps. This calculation 
required approximately 7 hours on a 25 MHz 486-based personal 
computer . 

Figures 5-7 show the magnitude of the Fourier transforms of the 
voltage and current results of Figures 2-4. The current results 
indicate the complicated frequency domain behavior of the coupled 
dipole system. 

The self admittance was obtained by dividing the complex Fourier 
transform of the driven dipole current by that of the Gaussian 
voltage pulse at each frequency. The results are shown in 
magnitude and phase in Figures 8 and 9 and compared with ESP4 
Moment Method results. Considering the differences in how the 
feed region is modeled (a 1 cm gap in the FDTD calculations vs an 
infinitesimal gap in ESP4) the agreement is quite good. 

The mutual admittance was obtained in a similar manner, dividing 
the complex passive dipole current by that of the Gaussian pulse. 
The results are shown in Figures 9 and 10. Again the agreement 
is quite good considering the different approximations and 
assumptions made in the FDTD approach relative to the ESP4 
computer code. 


IV Conclusions 

The capability of the FDTD method to predict mutual coupling 
between antennas was demonstrated. The test case was two 
parallel wire dipoles of different lengths, with one driven by a 
Gaussian pulse. The complex self and mutual admittance results 
obtained using FDTD showed good agreement with results obtained 
using the Method of Moments. 
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57 CM Dipole Feed Voltage 
spacing, 43 cm passive dipole, 61x 













57 CM Fed Dipole Current 
spacing, 43 cm passive dipole, 6 lx 












43 CM Passive Dipole Current 
spacing, 57 cm fed dipole, 61x51x80 Space 



Frequency 
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43 CM Passive Dipole 
spacing, 57 cm fed dipole, 61x51x80 Space 



Frequency 









43 CM Passive Dipole 

10.5 cm spacing, 57 cm fed dipole, 61x51x80 Space 








